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Abstract 

This paper deals with an optimal position management problem for a market maker 
who has to face uncertain customer order flows in an illiquid market, where the market 
maker’s continuous trading incurs a stochastic linear price impact. Although the exe¬ 
cution timing is uncertain, the market maker can also ask its OTC counterparties to 
transact a block trade without causing a direct price impact. We adopt quite generic 
stochastic processes of the securities, order flows, price impacts, quadratic penalties as 
well as security borrowing/lending rates. The solution of the market maker’s optimal 
position-management strategy is represented by a stochastic Hamilton-Jacobi-Bellman 
equation, which can be decomposed into three (one non-linear and two linear) backward 
stochastic differential equations (BSDEs). We provide the verification using the standard 
BSDE techniques for a single security case. Eor a multiple-security case, we make use of 
the connection of the non-linear BSDE to a special type of backward stochastic Riccati 
differential equation (BSRDE) whose properties were studied by Bismut (1976). We also 
propose a perturbative approximation scheme for the resultant BSRDE, which only re¬ 
quires a system of linear ODEs to be solved at each expansion order. Its justification and 
the convergence rate are also given. 

Keywords : BSDE, BSRDE, asymptotic expansion, portfolio, inventory, liquidity cost 
AMS subject classification: 91G80, 60H10, 93E20, 34E05 

1 Introduction 

The financial market currently being formed in the aftermath of the great financial crisis 
looks completely different from the previous one. Mandatory clearing for the standardized 
financial products and much higher regulatory costs for the rest of over-the-counter (OTC) 
contracts made many investors withdraw from the long-dated exotic derivative business and 
pay much attention to the trading of listed products with exchanges or standard contracts 
with central counterparties. 

In the new market, it is clear that the exchanges and central counterparties are the most 
important trading venues and have started to play a much bigger role than before. However, 
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these new developments have not completely diminished the importance of the traditional 
key players in the market, that is a market maker. A market maker is a firm that quotes buy 
and sell prices for financial securities and derivatives, and stands ready to perform these deals 
on a regular and/or continuous basis. They are crucially important to maintain liquidity for 
equities, currencies, commodities, government/corporate bonds, many structured products 
and derivatives. Even for products tradable at an exchange, market makers are playing 
an important role by intermediating non-financial corporates and other investors since it is 
not always possible for them to satisfy many regulatory conditions required to get a direct 
access to the exchange. There exist many other benefits such as those related to accounting, 
anonymity and flexibility that may be obtained with the help of market makers. 

Especially due to the proposed regulation on the leverage ratio and the higher capital 
amount required for the open positions, the market makers have to deal with formidable 
tasks. Due to the smaller warehousing capability of their balance sheets, they need more 
active position management. At the same time, they have to optimize execution strategies in 
order to avoid unnecessarily big market impacts and the associated transaction costs. See [8], 
for example, and other articles available in Risk.net to get some images of the recent market. 

In this paper, we consider the optimal position management problem for a market maker 
who is facing uncertain customer orders. We are interested in a good market maker who 
accepts every customer order with a predefined bid/offer spread. The spread can be stochastic 
but we do not allow the market maker to control its size dynamically based on its proprietary 
reasons in order to give a bias to the customer flows. Otherwise the firm will not be considered 
as a trustful market maker 0. We suppose that there exists a relatively liquid market for 
security borrowing and lending (i.e., so called repo transactions), which can be used by the 
market maker to answer the incoming customer orders. In addition to matching an incoming 
order directly to the security being warehoused in its balance sheet, it is assumed that the 
market maker can access two external trading venues. One is a traditional exchange where 
the market maker carries out absolutely continuous trading that incurs, however, a stochastic 
linear price impact. It is also supposed that the participants of the exchange can partially infer 
the inventory size of the market maker and that their aggregate reactions, as a preparation 
for the market maker’s future unwinding, affect proportionally to the security price. Another 
venue is the aggregate of the market maker’s OTC counterparties with which the hrm can 
execute a block trade without directly affecting the price in the exchange. In this case, 
however, the execution timing is uncertain. The modeling of the latter venue (we call it the 
dark pool) is closely related to the one introduced by Kratz & Schoneborn (2013) [3T] except 
that we allow stochasticity in its execution intensity. 

There now exist the vast literature on the optimal execution problems. Our model is 
closely related to the line of developments made by Bertsimas & Lo (1998) [9], Almgren 
& Chriss (1999, 2000) O 0], Schied & Schoneborn (2009) [1^ and to more recent works 
Ankirchner & Kruse (2013) [6], Ankirchner, Jeanblanc & Kruse (2014) [5] and Kratz & 
Schoneborn (2013) [3T]. In this first approach, the price process of the relevant security 
is exogenously modeled. There exist many other interesting approaches, such as models of 
supply curves ( See, for example. Bank &: Baum (2004) [7], Cetin, Jarrow &: Protter (2004) [TH] 
and Roch (2011) [39]. ) and those directly modeling the dynamics of Limit Order Books ( 
See, for example, Obizhaeva & Wang (2013) [36], Alfonsi, Truth & Schied (2008, 2010) [HIS], 
Truth & Schoneborn (2014) [20], and Cartea & Jaimungal (2015a,b) [HKTS]. ). We refer to 
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the review articles Gatheral & Schied (2013) [26] and Gokay, Roch & Soner (2011) [27] for 
the recent developments, various other aspects and references. 

The two main differences of the current work are the focus on the market maker’s posi¬ 
tion management problem with uncertain customer orders, and the generality in its setup. 
We allow quite generic stochastic processes for the securities, position impact factors, com¬ 
pensators of incoming customer orders, execution intensities of the dark pools, repo rates 
relevant for the security borrowing/lending, and the quadratic penalties for the outstanding 
position size etc. To the best of our knowledge, it is the most generic setup among the lit¬ 
erature adopting the first approach in the last paragraph. The resultant optimal strategy 
becomes fully adapted to the market filtration instead of a deterministic one usually found 
in the literature. In contrast to a very short term liquidation strategy for which non-random 
(or even constant) coefficients may be sufficient, this generality is necessary for the position 
management problem with a longer time horizon in which a significant change of the mar¬ 
ket conditions is naturally expected. Although the first approach we choose is a simplified 
reduced-form approximation of Limit Order Books, it allows more flexible modeling of the 
underlying processes including multiple securities and also their mutual dependence, which 
is expected to be more relevant for our medium term problem. 

We follow the technique proposed by Mania & Tevzadze (2003) [35] to derive the relevant 
stochastic HJB equation and its decomposition into three (one non-linear and two linear) 
backward stochastic differential equations (BSDEs). For a single security case, we use the 
standard results on BSDEs and the comparison theorem (See, for example, Ma & Yong 
(2007) [33] and Pardoux & Rascanu (2014) [57].) for verifying the solution. For a multiple- 
security case, however, we need to handle a matrix-valued non-linear BSDE for which we do 
not have an appropriate comparison theorem. We show that the relevant BSDE is actually a 
special type of backward stochastic Riccati differential equations (BSRDEs) associated with a 
stochastic linear quadratic control (SLQC) problem. Interestingly, a seemingly quite different 
setup of the optimization problem gives rise to the same BSDE. Thanks to this relation, we 
can guarantee the existence of a uniformly bounded solution by theorems proved by Bismut 
(1976) [12]. The main difficulty for the implementation of the proposed scheme is the concrete 
evaluation of this BSRDE. We propose a perturbative expansion technique for the BSRDE 
with a general Markovian factor process, which only requires to solve a system of linear ODEs 
at each order of expansion. A justification and convergence rate of the approximation scheme 
are also given. 

The organization of the paper is as follows: Sections 2 and 3 give some preliminaries, the 
detailed market description and the market maker’s problem. Section 4 gives the derivation of 
the candidate solution and its verification. An extension to a multiple-security case is given in 
Sections 5 and 6. Sections 7 and 8 deal with implementation. In particular, the perturbative 
scheme and its error estimate are given in Section 8. The behavior of the terminal position 
size with respect to the penalty size is studied in Appendix A. 

2 Preliminaries 

We consider a complete filtered probability space, in which all the stochastic processes are 
defined, (D,7^, F, P) where F = (JG)i>o is the filtration satisfying the usual conditions. W is 
the d-dimensional standard Brownian motion and the P-augmented filtration generated by 
W is denoted by F'^ = )t>o- We assume that F’^ is a subset of the full filtration; F^ C F. 
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For the ease of discussion, let us define the following spaces of the stochastic processes (p > 1): 
• E>rit,T) is the set of progressively measurable process X taking values in and satisfying 


E 


i^iirt,T] 


:= E 


sup 

■sG[t,T] 


< oo 


where we use the notation 

lkll[a,fe] := sup{|xt|,t G [a,b]} 

for X : [0,T] ^ M’’. We write ||x||[o^t] = ||a;||t- Its norm is defined by 

i/p 


X 


:= {E[||X||f,_^J} 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


Mr{t,T) is the set of progressively measurable process X taking values in M’’ and satisfying 

(2.4) 


E 


rT \ P/ 2 ' 

iXA^dt] 


J 


< oo, 


and its norm is defined by 




\XA^dsj 


p/2 'i l/P 


(2.5) 


In every space, the subscript r may be omitted if the associated dimension is clearly seen 
from the context. 


3 A single security case 

Firstly, let us summarize the standing assumptions. They are obviously not the weakest ones 
but allow simple analysis and also do not make the model unrealistic in a practical setup. 
Note that the definition of each variable will appear along the discussions in the following 
sections. 

Assumption A 

M{uj, dt, dz) is a random counting measure of a marked point proeess with a bounded support 
K C E\{0} for its mark z, and H is a counting process. All the stoehastic processes whieh 
do not jump by Af and H are assumed to be -adapted and hence continuous. This E^ 
adaptedness ineludes all the stochastic processes defined below. 

(oi) 5 ; n X [0, T] —)■ M is non-negative and S G §^(0,T). 

(02) 6, f : n X [0, T] —)■ R and b, I G S^(0, T). 

(03) A(-, •) : n X [0, T] X R —)■ R is such that A(t, •)(w) is a non-negative measurable funetion 
with bounded support K C R\{0} for every t G [0,T] and w G If, and that A(-,z) is a uni¬ 
formly bounded -adapted process for every z £ K. 

(04) 7 : n X [0,T] ^ R are uniformly bounded and non-negative. 


4 



(as) M,rj,X ; O x [0, T] —)■ M are uniformly bounded and strictly positive. 
(ae) ^ : O —)■ M is strictly positive, bounded and -measurable. 

(ay) /3 : X [0, T] —)■ M is uniformly bounded. 

(as) There is no simultaneous jump between M and H. 

Notation : For a bounded variable x, we denote its upper bound by x. 


3.1 The market description 

We are interested in a market maker, who has to face uncertain customer orders regarding the 
single specified security. An extension to a portfolio management including multiple securities 
will be discussed in later sections. As a good market maker, the firm accepts every customer 
order with a predefined bid-offer spread. Although the spread can be dynamic depending on 
the external market variables such as the security’s volatility, it is supposed that the firm does 
not adjust it in order to control the customer flows based on the firm’s proprietary reasons. 

The market maker is assumed to buy and sell the security through the two major trading 
venues. The first venue is a standard exchange, where the market maker carries out absolutely 
continuous trading. It is assumed, however, to incur a linear stochastic price impact. In 
addition, the participants of the exchange (partially) infer the inventory size of the market 
maker. They expect future buy/sell orders from it and adjust their positioning accordingly. 
We assume that the aggregate effects of the participants change the market price by the 
amount proportional to the inventory size, which becomes another source of the market 
maker’s future trading costs. 

The second venue is the aggregate of the OTC block trades with the firm’s customers or 
the dark pools. The market maker can buy/sell a block trade without directly affecting the 
market price, however, its timing is assumed to be uncertain. Although we call this venue the 
dark pool, it actually means the aggregate of OTC block trades with the firm’s counterparties 
as well as potentially multiple dark pools to which the market maker can access. It is 
a simplistic model for which we do not consider an order-size dependent intensity process 
nor possibility of the partial execution. Unfortunately, this seems unavoidable to keep the 
problem tractable. 

In addition to the above two trading venues, the market maker can match an incoming 
customer order to its outstanding position being warehoused in its balance sheet. This is the 
distinguishing feature of the market maker. Because this is the most profitable way to reduce 
the position, the market maker needs to adjust buy/sell orders based on the expected future 
customer flows. If the market maker cannot answer an incoming customer order within its 
inventory, it needs to borrow the security through the corresponding repo market by paying 
the stochastic repo rate. On the other hand, when its inventory is positive, the market maker 
earns money by lending the security through the repo market. 


We model the the market-maker’s position at time s > t starting from the position size 
X G M at time t as 


Xf’^{t, x) = X -\- 



f'S ps 

zM{du, dz) + / TTudu / 5udHu 

Jt Jt 


(3.1) 


It JK 

where the second term describes the customer flow, which is represented by the marked point 
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process expressed by the counting measure J\f. K C K\{0} is a bounded support for the mark 
z which gives the size and direction (+ or —) of the order d. (vr, (i) denotes an F-predictable 
trading strategy of the market maker through the exchange and the dark pool, respectively. 
H is the counting process, whose jump signals the happening of an execution event in the 
dark pool. For simplicity, we assume no simultaneous jump between J\f and H. The negative 
position size X'^’^ < 0 is always interpreted as a short position taken by the security borrowing 
through the repo market. Let us assume the existence of the compensators A for M (A for 
H) so that 


J\f {ds, dz) = j J (^Af{ds,dz) — A{s, z)dzds^ 


0 JK 

dH,^ = I ( dHg — Xsds 


/ 


(3.2) 


for t G [0, T] are F-martingales. This also implies that an occurrence of a customer order and 
an execution in the dark pool are totally inaccessible. For later convenience, let us define 




<l>i := / zA{t,z)dz, 'kt := / \z\A{t,z)dz, ^ 2 ,t ■= / z‘^A{t,z)dz (3.3) 


IK 


L 


for t G [0,r], which are the moments of the size of the customer orders. By Assumption A, 
the above processes are uniformly bounded. 

The price observed in the exchange S’^’'^(t,x) i.e., the market price under the impact of 
the market maker’s strategy (tt, 5) starting from the position size x at time t, is assumed to 
be given by 

x) = Ss + MsTTs - x) (3.4) 


for s G [t,T]. The second term denotes the stochastic linear price impact, where M is the 
F^-adapted impact factor. The last term denotes the aggregate impact from the market 
participants’ reactions to the market maker’s inventory size. 

Notice that we are not assuming the perfect observability of the market maker’s position 
X to the other investors. It is likely that they can infer X only vaguely. Thus, their reac¬ 
tions likely have a big noise, too. However, this noise part can easily be absorbed into the 
definition of S', the unaffected price of the security. For the market maker’s point of view, 
X is directly observable and /3 is simply its coefficient which can be obtained by the linear 
regression of the security price. We model (J as a uniformly bounded process possibly being 
correlated with other market variables, such as volatility of the security. Due to the presence 
of the customer orders, the last term is not directly determined by the trading volume in 
the past and hence different from the standard model of the permanent price impact. It is 
more closely related to the works on the large trader’s problem studied by Jarrow (1992) [29], 
Cvitanic Sz Ma (1996) |18] . and Bank & Baum (2004) [7]. 


We model the cash flow in the interval ]t,T] to the market maker with strategy (tt,;!) in 
^This boundedness can be interpreted as the maximum acceptable order size set by the market maker. 
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the following way: 


f Sg’^{t,x)TTsds — f [ S^!^{t,x){l —sg'a.{z)bs)zM{ds^dz) 

Jt Jt Jk 

-j ({Ss-l3sX::!{t,x))5s + %%?)dHs+ j lsX^/{t,x)ds . (3.5) 


Let us explain the economic meaning of each term below: 

• (1st term) The cash flow from the trades through the exchange. 

• (2nd term) The cash flow from accepting the customer orders with a (proportional) bid/offer 
spread b. 

• (3rd term) The cash flow from the trades through the dark pool. 

• (4th term) The cash flow from the security borrowing/lending with a repo rate 1. 


We need additional comments for the third term describing the trades with the dark pool. 
Firstly, the basic transaction price is given by 

Ss-psX::_\t,x) (3.6) 

which does not include the price impact from the continuous trading of the market maker. 
Inclusion of MgTTs to the price could induce price manipulation, and more importantly, the 
trading counterparties will not accept expensive price caused by the market maker’s temporal 
trading activity. We also add the spread p\5\ to the above priced as a premium that the 
market maker pays to the counterparty who has accepted a block trade. 

We consider T^l (year) as a relevant time span for the control of the market maker. 
More realistically, it can be a Quarter or a half year, and we neglect the net proceeds from 
a money market account for this time interval. This can be understood as a (nearly) zero 
interest rate, or equivalently, we can interpret that the cost function (see below) is given in 
the discounted basis. We are also interested in relatively liquid market in a sense that the 
borrowing and lending of the security is always possible as long as a given stochastic repo 
rate is paid. In a highly illiquid market, neither seamless execution of market orders in the 
exchange nor a functioning repo market can be expected. 


3.2 The market maker’s problem 

Definition 3.1. 0 We define the admissible strategies U by the set of ¥-predictable processes 
(7r,(5) that belong to BI^(0,T) x BI^(0, T) and also Markovian with respect to the position size, 
i.e., they are expressed with some measurable functions {fi^, f^) by 

TT, = r{s,x:f{t,x)), 5 , = fis,x::_\t,x)) (3.7) 

where, for a € {vr, (i}, /“ : x [0, T] x M —)• M and f^{-,x) is an ¥^-adapted process for all 
x gR. 

®Thus the additional cost is given by rj\3\^. 

^ It may not necessary to constrain the admissible strategies as Markovian with respect to X. However, 
limiting the strategy space at this stage makes the following analysis much clearer. 
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We suppose that the market maker tries to solve the following optimization problem: 


V(t, x) = ess inf E 
{n,5)£U 


+ [ {Ss'^{t,x)Trs - lsX'^'^{t,x))ds + [ [ Sg!^{t,x){l - sgn{z)bs)zN'{ds,dz) 

Jt Jt JK 

+ r ([Ss - /3 sX:!_^( t, x)] ,5, + 15,12) dH, 


(3.8) 


The first two terms are introduced to give penalties for the outstanding position size. It 
is natural to consider that ^ and 7 are proportional to the variance of the price process of 
the security. One may also want to take into account the regulatory costs arising from the 
outstanding position in the balance sheet. It is possible by an appropriate modification of 7 
and I as long as the relevant costs can be reasonably approximated by a quadratic function 
with respect to the position size X. Note that the coefficients of the quadratic function can 
be stochastic. 


We can observe that the expectation in (13.811 is finite for all ( 7 r,( 5 ) G U. This can be 
easily checked by the fact X'^'^{t,x) G §^(t,T) and S^’^{t,x) G However, due to the 

2 nd order terms of {7r,6) arising from {—(dX'^'^Tr) and {—(dX'^'^S), the cost function could be 
unbounded from below, and then the problem would be ill-defined. In order to guarantee the 
well-posedness of the problem, we need additional assumptions. 

Firstly, let us write the dynamics of the F^-adapted bounded process /3 as 

dfit = Hfdt + a^dWt ■ (3.9) 


Furthermore, we denote 


e:=e- 


h 

2 ’ 



~ I ^ 

'i-'t + Y- 


(3.10) 


Assumption B 

( 61 ) : n X [0, T] —>■ M, X [0,T] are uniformly bounded, and -adapted. 

( 62 ) 7 is non-negative dF ® dt-a.e.. 

( 63 ) There exists a constant c > 0 such that > c a.s. and M,r] > c dP® dt-a.e.. 

Definition 3.2. The cost function for the market maker with a given position size x G M at 
t G [0, T] is 


J‘’*(7r,d) =E 

rT 


f\X^/{t,x)\^+ f (^^s\Xf^\t,x)\^ + Xf^\t,x){Psbs^,-ls))ds 


+ / + ^sVs^s + + MsQs)t^s + iSsAiids] -|- [SgQg -|- ^^2,s)'jds 


J-W 

•'t 


(3.11) 


where 0 : If x [ 0 , T] —>• M is defined by 0^ := — bg^g. 
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Proposition 3.1. Under Assumptions A and B, the market maker’s problem h3.^) is equiv¬ 
alent to 


P(t, x) = ess inf (3-12) 

('7r,5)EZV 

and it has a unique optimal solution {7r*,6*) G U. 

Proof. Applying Ito-formnla, one obtains 

It ^ ^ I Y^sdHs 

+ f f Y^‘^d\f{ds,dz)+ f ^\Xf’\t,x)\'^{n^,ds + cr^dWs) . 

J t- J K fj t 

Then, replacing the /^-proportional terms in (|3.8I) fdsXf:^{t, x)[7rsds + 6sdHs]^ by using 

the above relation and p.lOl) yields 


V(t,x) — ^x"^ = ess inf E 

2 {■K,&)elA 


i\X^'\t,x)\'^ + f (js\x:’\t,x)\^ + Xf’^iPxXPsbs'^s - ls))ds 


+ J ■|Afs 7 r^ + + [(< 5 s + MsQs)TTg -)- + {Sg&s + ”^^2,8) ds 


Tt 


(3.13) 


where the integrals by the counting measures are replaced by their compensators. Here, 
one can check that the local martingales are true martingales under the assumptions. In 
particular, one can use the Burkholder-Davis-Gundy (BDG) inequality, the fact that G 
§f{t,T) and the boundedness of ap for the dW integration term. For the jump part, it 
suffices to check that the integration by the corresponding compensator is in L^(H) (See, for 
example, Corollary C4, Chapter VIII in [l3].), which can be confirmed by the boundedness 
of the compensators, A and A. 

Because all the processes except {X’’^’^, vr, 6) are F'^-adapted and (vr, 5) satisfies (13.7p . 

the expectation conditioned on T) in (j3.13p can be replaced by V thanks to the 

Markovian nature of . Notice that the information of the counting measures M, H only 
appears through the position size X’’^’^. However, Xf’^{t,x) = x has been already fixed. 

Thus, redefining the value function by V{t, x) := V{t, x) —obtains the result (|3.12p 
as an equivalent problem for the maker maker. 

The remaining claims easily follow from the standard arguments (See, for example, The¬ 
orem 3.1 in Bismut (1976) [E].) since now all the quadratic terms have positive coefficients. 
For simplicity, let us consider the case where the initial time is zero, t = 0. The cost function 
is a continuous map from to M and obviously strictly convex. It is also proper since, 

for example, J°’^(0) < oo. We also have the so-called coerciveness since 

J*^’^(n)oo, when | |ti| oo . (3.14) 

The above observations and the fact that U is a Hilbert space tell us that, for a large enough 
a G M, the set {rt G (u) < a} is non-empty, convex and weakly-compact. Thus, there 
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exists an minimizer, which is unique due to the strict convexity of the cost function. 


□ 


Remarks on in ()3.4p 

The presence of term in (j3.4l) is not necessarily appropriate for every type of investors. 

For example, suppose that the investor is risk-neutral and /3 is positive. In this case, the 
investor may accumulate an extremely large long position which would make the security 
price significantly negative. The investor can receive positive cash flow by further increasing 
her long position which makes the system ill-defined. However, as we have seen in the above 
discussion, it does not cause any regularity problem under mild conditions regarding the 
penalty size on the outstanding position of the investor. Although one may feel uneasy by 
the fact that the well-posedness of the model depends on the risk-averseness of the agent, 
we think that this term makes the model more realistic for the market maker. In fact, this 
term is expected to arise exactly because the other investors know that the relevant market 
maker has to operate with a rather stringent position limit i.e., risk averse. From the view 
point of the market maker, it is being squeezed by the other investors as long as there exists 
an information leak about its position size. 


4 Solving the problem 

4.1 A candidate solution 


Let us prepare the optimality principle for the current problem. 


Proposition 4.1. (Optimality Principle) Let Assumptions A and B are satisfied. Then, 
(a) For all x G M, (vr, 6) gU and t G [0, T], the process 


^ V{s,X^’^{t,x)) + (ju\X((’\t,x)\^ + X0’\t,x){l3ubu'^u-lu)^du 

+ J + ^uVud'^ + [(S'm + MuQu)'^u + Su^udu] + {Su®u + 

is an ¥-submartingale. 

(b) is optimal if and only if 


se[t,r] 


^ V{s,Xf’^*{t,x)) + {t,x)\‘^ + Xf’^* {t,x){/3ubu^u - lu)^du 

+ J (j^uT^*u T ^u'Hud*u + [{Su + u + Su)^ub*u\ + {Su&u H- ^^2,u)^du 

is an ¥-martingale. 


se[t,T] 


Proof. One can easily confirm it from the definition of the value function V, the fact that 
V{T,X^’^) = and the form of the cost function See, for example, Propo¬ 
sition (A.l) of Mania &: Tevzadze (2003) [35]. □ 
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Firstly, by following the method proposed by Mania &: Tevzadze [35], we derive the BSDEs 

from the necessary condition so that the above optimality principle is satisfied. Then, we 

are going to show that there exists a solution for every BSDE and confirm that it actually 

satisfies the optimality principle. This gives us one optimal solution. But we know that the 

solution is also unique due to Proposition 13.11 

Let us assume that the semimartingale ( V(t,x )) has the following decomposi- 

V /te[o,r] 

tion for every x G M: 


V{s,x) 


rs ps 

V{t^x)+ / a{u^x)du+ / Z{u^x)dWu 

Jt Jt 


(4.1) 


where a : x [0, T] x M ^ M, Z : H x [0,T] x M ^ and a(-,x) as well as Z(-,x) are 
F^-adapted processes for all x G M. Let us suppose V{t,x) are twice differentiable with 
respect to x. By applying Ito-Ventzell formula, we obtain 


PS PS 

V{s^X'^^^{t^x)) = V{t^x) + / a{u^X^^^{t^x))du / Z{u^X^^^{t^x))dW^ 

Jt Jt 

+ V^{u,X^'^{t,x))7rudu +j^(v{u,X^’^{t,x) + z) 

+ (Lx) + <5„) - V{u,XZ:^it,x))yHu . 


'u, dz) 
(4.2) 


Separating the local martingale parts, a necessary condition for the optimality principle is 
given by 


a{u,x) + J \^V{u,x + z)-V{u,x)jA{u,z)dz+ juX + xiPubu^u-D 
+ (*S'„0„ + + inf < 14 (w,x) 7 r + {V{u,x + 6) - V{u,x))Xu 


+iVLi7r^ + + {Su + MuQu)Tr + SuXud z — 0 


(4.3) 


dP 0 dt-a.e. in fl x [0, T] for every x G M. 

Substituting the resultant drift term a(-, •) into (|4.1I1 yields a backward stochastic PDE 
V{t, x) = l^j[V{u,x + z)-V{u, x)] A{u, z)dz + 7 „x^ + x{Pubu^u - L) 

+ (*S'„0„ + I + [yiu,x + 6) - V{u,x)]Xu 

+MuTt‘^ + XurjuS"^ + {Su + MuQu)t^ + 5'«A„d| ^ x)dWu, (4.4) 

which is sometimes called a stochastic HJB equation, (vr, <5) should be chosen for each u G 
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[t, T], Exploiting the quadratic nature, let us hypothesize that, for every t G [0, T] and x G M, 


V{t,x) = V 2 {t)x^ + 2Vi{t)x + Vo{t) 
Z{t, x) = Z 2 {t)x^ + 2Zi{t)x + Zo(t) 


(4.5) 

(4.6) 


where V 2 , Vi, Vq ^ x [0, E] ^ M and Z 2 , Zj, Zo : Q x [0,T] are F'^-adapted processes. 

Then, (14.3j) can be rewritten as 

a{u,x) + { 2 V 2 {u)^uX + V 2 {u)^ 2 ,u + 2Ei(u)4>„) + + x(/3„6„T„ - /„) + (^Su&u + ^^ 2 , 

[V 2 {u)x + Vi{u) + ^{Su + MuQu)] ^2 

M„, 


+ inf < Mu (vr + 
n,5 I 


[^(^) + Vu] + 


[^ 2 (^) 3 : + Vl{u) + 

E2(tt) + rju 


= 0 d¥ ® dt — a.e.. 

For the well-posedness, we must have V 2 + r] > 0 dF ® dt-a.e.. 

Gathering each of (x^, x^, x*^)-proportional terms in (14.4p . one obtains the following result. 


(4 


A Candidate Solution 

A “candidate” of the optimal solution and the corresponding value function for the market 
maker’s problem hS.12\) are given by 


K = -^(v 2 iu)x^_’^ {t,x) + Vi{u) + ^{Su + Mue^ 


1 


< 5 : = - 


[V2iu)x:_’^ (t,x) + Vi{u) + -Su] 

V2{u) + r]u 


(4.8) 


(4.9) 


for u G [t,T] and V{t,x) = V 2 {t)x‘^ + 2Vi{t)x + Vo{t), respectively. Here, X'’’'*(t, x) is the 
solution of 


Xf (t, x) = X + / / z. 

It Jk 


ns fS 

'M{du,dz) + / Tr*du+ / 6 *dHu, sG [t,T] . (4.10) 

Jt Jt 


(V 2 , Z 2 ), (El, Zi) and {Vq, Zq) must be the well-defined solutions of the following three BSDEs 


•i-(4- + 


A 7 , 


t { \Mu V2{u)+r]u 


)V2{uf +lu^du- Z2{u)dWu (4.11) 
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VM = -l + 


+ V2iu) 


1 


+ 


A 7 , 


iMu V2{u)+r]ui 2 2 


— - ^ 0 n - 


du- J Zi{u)dWu (4.12) 




■«) + ^) - li(it)(4>n + 


-V2{u)^2,u - + Pu^2,u) + ^MuOl ^du- Zo{u)dWu , (4.13) 

satisfying 

dP <Si dt-a.e. mflx[0,T]. 


V2 + r] > 0 


(4.14) 


4.2 Verification 


We are now going to study each BSDE and show the existence of the candidate solution, and 
also confirm that it actually satisfies the optimality principle. 


Proposition 4.2. Under Assumptions A and B, the BSDE |.^.11[ ) has a unique solution in 
{V 2 , Z 2 ) € S^(0, T) X BI^(0, T) for V p > 1, and in particular ¥ 2 ( 1 ) satisfies for every t G [0, T] 
and e > 0 that: 


E 

(7 + 

r(w 

+ —'jds 




It 

hs' 



< 


(T-t + e) 


rE 




T < Mt) 

T 

J (yAIg + (T — s + e)^ 7 s ^ds 




w 


(4.15) 


Proof. Let us define the function as 

/(i>y) =-(^ + + 7t ■ (4.16) 

\Mt y + rit^ 

Firstly, let us consider the BSDE 

f{s,Y,yD)ds- ZsdWs (4.17) 

Due to Assumptions A, B and the definitions (I3.10p . ^ and f{t,y V 0) with any fixed y are 
bounded. Furthermore, it is clear that f{t, y V 0) is a decreasing function in y. Thu^ (|4.17l) 
satisfies the standard monotone conditions for the BSDE. By Theorem 5.27 in m □, there 
exists a unique solution {Y,Z) G §^(0,T) x ]Hl2(0,T) for all p > 1. On the other hand, it is 

®One can simply put /i—l—O in the theorem. 
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clear that we have a trivial solution (Y, Z) = (0, 0) if .^ = 0 and 7 = 0. Since the terminal 
value and the driver / is increasing in (■^, 7 ), we actually have > 0 by the comparison 
theorem (See, for example, Proposition 5.33 in [37].). As a result, the BSDE (14.111) has a 
unique solution {¥ 2 -, Z 2 ) G S^(0, T) xBI^(0, T) for all p > 1, and in addition, V 2 is non-negative. 

The derivation of the upper and lower bounds is an adaptation of Proposition 2.1 in 
Ankirchner, Jeanblanc &: Kruse (2014) [5] for our problem. Let us start from the derivation 
of the upper bound. For all y, G M, 


— 2ky + k'^ > 0 


(4.18) 


is satisfied. For an arbitrary constant e > 0, choosing k = 


M, 




\Mt y + rjt 


2 1 2 2 Mt 

^ - “r-f + e^^ (r-t + e )2 


T-t + e 

y + 


yields 


(4.19) 


for all y > 0 . 

With some abuse of notation, consider the next linear BSDF 


Yt^ = i + 


j: 


T-s + e 


y: + 


M, 


(T-s + e) 


+ 7 s ^ ds - 


ZtdW, 


(4.20) 


This is a linear BSDE with a bounded Lipschitz constant. Due to the boundedness of M, 7 , 
there exists a unique solution (W,Z^) G §P{0,T) x ]Hl2(0,T) for all p > 1. By the inequality 
(I4.19P and the comparison theorem, we have 


V2{t) < Yl 

for all t G [0,T] and e > 0. In addition, Y^ can be solved as 


(4.21) 


Yf = E 


Ce 






w 


(T-t + e)2 


E 


e^^ + J (^Mg + (T — s + e)^7s^ds 


-jrW 


(4.22) 


and hence we obtained the desired upper bound. 
Now, let us study the lower bound. Put 


Ft :=E 


yI 


^ ' I As 


( 


VMs rjs/ 


ds 




w 


Due to the existence of a constant c > 0 such that M, y > c, it satisfies 

-= Y V <—^1 + T(1 + A)^ [:= k ) 


(4.23) 


(4.24) 


where A are the upper bounds of A, respectively. Therefore, there exists Z G Hj(0, T), Vp > 
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0 such that 


dVt = - 


—+ ^ 
Mt rjt 


dt + ZfdWt . 


Then the process Ut := 1/Vt, which has the bounds 1/k <U < ^, satisfies 


Ut = ^ + 


J_ ^ 
Ms rjsJ 




Us 




TsdWs , 


(4.25) 


(4.26) 


where T := — Here, the terminal value, the first term of the driver and the coefficient of 

|r|2 are all bounded. Thus the comparison theorem for the quadratic BSDE (See, Theorem 
2.6 in [30]), one sees Ut < V 2 {t) and hence the desired result is obtained. □ 

Since the BSDEs for (Vi,Zi) and (Vo,Zq) are linear, one can use popular established 
results to obtain the next Proposition. 


Proposition 4.3. Under Assumptions A and B, there exist unique solutions (Vi,Zi) E 
S^(0,T) X 11^(0,T) for and {Vq,Zq) E S^(0, T) x ]HI^(0,T) for respectively. 

Proof. We denote by C some positive constant, which may change line by line. Erom Propo¬ 
sition [421 V 2 is uniformly bounded and hence so is the linear coefficient of Vi. In addition. 


E 


\{f 

V2{S)(' 


\ L 


Ms 


+ 


Xs 


V2{s) + rjsi 2 


5, 0 . 


< CE 


1 + / ( )<^'S 


< 00 


\ 1 . 


bs'its) - -{fdsbs'^s - h) 

ds j 


(4.27) 


Thus, by Theorem 5.21 (see also Section 5.3.5) in [37|, there exists a unique solution (Pi, Zi) E 

§^(0,r) X ]H4(0,T). 

As for ()4.13p . it is easy to see that 


E 


1 


+ 


A. 


Ms V2{s) + r]s 


—-('S's0s + Ps^2,s) + -^MsQI 


S) + y)' - Vt{s){^s + bs^s) - V2(s)^2„ 


ds 


< CE 


1 +y ^ihi('S)i^+issi^+ 


< 00 , 


(4.28) 


where we have used Pi E §'^(0,T) proved in the previous arguments. Thus, by the same 
reasoning, there exists a unique solution {Vq,Zq) E §^(0,T) x ]HI^(0,T). □ 


In order to check the optimality condition, we also need the following property of X'^ 


Proposition 4.4. Under Assumptions A and B, the process of the position size 
given by belongs to §^(t,T). 



se[t,r] 
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TT* 5* 

Proof. Let us take the starting time 0 and write Xg ’ (0, x) as Xf for simplicity. Then, 


Xg ’ {0,x) = X + I I z. 

'o Jk 


{V2{u)x:y\o,x)+viiu) + f) 


■M{du, dz) — [ '' —-— \ \ " — "" — —dH„ 

Jo V2[U) + Vu 

r I + V(\\ )xf’'\o,x) + (^ + TTTT^) 

Jo + VM„ V2(u) + riu/\ 


\du 


(4.29) 


Under Assumptions A and B, there exists some positive constant C such that 


\x;\^ < c 


+ 


1 + / (iv;!* + |Vi{!,)|-‘ + ISJ* + \bp)du 

^Jo Jk ^ \ Jo V2{u) + r]u 


dHu 


for every t G [0, T]. Let us define a sequence of F-stopping times (r„)„>o by 

Tn ■= inf {t > 0 : 1X^*1 > n} A T , 


(4.30) 


(4.31) 


and denote the r^-stopped process of the position size as Since we already 

know that X* G S^(0, T), it is clear that Tn ^ T a.s. as n —>■ oo. 

The BDG inequality (see, for example. Theorem 10.36 in |28] for general local martingales) 
and positivity of integrand yield 


E 


IX 


*Tn I 4 


< CE 


1 + + |Ui(u)|4 + + \buf)du 


+ \ I J z^M{du,dz)j “*" ( y 


V2iu)X*J- + Ui(u) + % 


V2{u) + r]u 


dHn 


(4.32) 


Again by the BDG inequality and the boundedness of A, one has 


E 


* \ 2 
\x:j^\‘^dH, 


< 2E 


^ \ 2 
\x:j^\‘^dH, 


< CE 

One obtains, by similar analysis, that 


+ 2E 


(£|X-"|2A„dn)' 


[ \x*jydu 

Uo 


(4.33) 


E 


V 2 iu)X*'[l^ + Ui(n) + ^ 


y2{u) + r]u 


dHn 


< GE 


J^(\Vi{uy + \Su\^)du + \X*J-\^du 


(4.34) 
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Thus, it can be shown from (14.321) and the boundedness of K that 


|4 

1 t \ 

< CE 




1 + II^iIIt + II'S'IIt + 


E 


and hence, by the Gronwall lemma, for Vt G [0,T], 


\t + 


f 


\x:r^\^du 


E 


I I 4 


< CE [1 + \\Vi\\^ + \\S\\^ + \\b\\^] 


y] < oo . 


(4.35) 


(4.36) 


Passing to the limit n ^ oo, we see E 


\xt\^ 


< C for every t G [0,T] with some positive 


constant C. Using the BDG inequality and the above estimate, we obtain from (|4.30l) that 

pT 


E 


I 114 

I 11 T 


< CE 


1 + ll^illr + 11‘S'llr + ll^llr + / 

Jo 


< oo . 


(4.37) 

□ 


Corollary 4.1. Under Assumptions A and B, the candidate solution {7r*,6*) given by 
and \4.S^ is well-defined, unique and satisfies {Tr*,S*) G S^(t, T) x S‘^(t,T) C U. 

Finally, we arrived the first main result of the paper. 

Theorem 4.1. Under Assumptions A and B, the candidate solution given by 

and is, in fact, the unique optimal solution of the market maker’s problem given by 

(KM- 

Proof. It suffices to confirm that the optimality principle of ProDosition l4.1l is indeed satisfied. 
Firstly, we have to see 


/ Ziu,Xf^^*{t,x))dWu 

't / sG[i,T] 

/7u(n, x^y* it, x)+6:) - v(u, x^y* (t, x))] dSu 

't ^ ^ / se[t,T] 

f f iviu,Xiy it,x) + z)-Viu,Xiy it,x))\Midu,dz) 

h Jk^ ^ / se[t,T] 


are all true F-martingales. For notational simplicity, let us put t = 0 and X* = Xf ’ (0,x). 

By the BDG inequality. Proposition 14.2114.31 and 14.41 there exists a positive constant C 
such that 




l-S 


/ rT \ 2 

E 

sup 

/ Ziu,Xf)dWu 

< CE 

/ \zis,x:yds) 


_sg[0,T] 

Jo 


\Jo ) 


< CE 


< CE 


pT \h / rT 

I 2 I \r^ I 2 , 


1^ \z2isy\x:\^dsj \ziisy\x:\^dsj \zoisyds 


l + ||X*||^ 


+ J ^ 1 ^ 2 ( 5 )^ + |Zi(s)|^ + |Zo(s)|^^ds 


< 00 


(4.38) 
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where in the second inequality, we have used the fact that 


(fl + 6 + ^ -|- \/b -|- y/c . 


(4.39) 


for every a,b,c > 0. Thus, ^ Z{u,X*)dWu^ ^ ^ is a martingale. 


For the integrations by the counting and marked point processes, it is suffice to check that 
the integration by the corresponding compensator is in L^(n) (See, Corollary C4, Chapter 
VIII in [l3].). Therefore, for the second term, we need to check 


E 


UO 


\Viu,X* + 6*J-Viu,X*)\Xudu <oo 


In fact, 


E 

L^o 
= E 


\viu,x:+6:)-v{u,x:)\x^du 


LJO 


Xu\v2iu){2x:s: + ( 5 :) 2 )+ 2v,{u)6:\du 


< CE 


L4o 


{\x:\‘^ + \6:\^ + \v,{u)\‘^)du 


< oo. 


Similarly, for the third term, 

’^[f I 

= E 


< CE 


\Viu,X* + z)-V{u,X*)\A{u,z)dudz 

[ f \V2 {u){2X*z + z"^) + 2 Vi{u)z\A{u, z)dudz 

Jo Jk 

£{\x:\^ + \V,{u)\‘^ + \<^2,u\)du 


< oo 


(4.40) 


(4.41) 


(4.42) 


where we have used the boundedness of the compensator and the support K. The above 
facts combined with the construction of a{t,x), strict positivity of M as well as X[V 2 + rj], 
guarantee that the optimality principle in Proposition 14.11 is indeed satisfied. □ 


In Appendix A, an investigation of the relation between the penalty size and the remaining 
position at the terminal time is given. It is proved that the terminal position size X^ can 
be made arbitrary small by increasing the size of the penalty This result implies that the 
proposed strategy can also be used for the liquidation problem in the presence of uncertain 
customer order flows. 


5 An extension to a portfolio position management 

In the following sections, we are going to extend the previous framework so that we can deal 
with the optimal position management for a market maker in the presence of n € N securities. 
Firstly, Let us summarize the assumptions below. As before, the definition of each variable 
will appear along the discussions in the following sections. 
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Assumption A' 

M^{uj,dt,dz),i G {1, • • • ,n} is a random counting measure of a marked point process with a 
bounded support K C K\{0} for its mark z, and H^,i G {1, • • • , n} is a counting process. All 
the stochastic processes which do not jump by AA*, H'^ for i G {1, •'' assumed to be 

-adapted and hence continuous. This adaptedness includes all the stochastic processes 
defined below. 

{a'fi) S : O X [0,T] —>• is non-negative and S G S^(0,T). 

(4) 6, Z : r? X [0, T] ^ and b, I G §^(0, T). 

K) X [0, T] X M — 7 > R for i G {1, • • • , n} are such that A*(Z, ■){io) is a non-negative 

measurable function with bounded support K C R\{0} for every t G [0, T] and w G fZ, and 
that A*(-, 2 ;) is a uniformly bounded -adapted process for every z ^ K. 

{a'fi) 7 : O X [0, T] —)• is uniformly bounded and takes values in the space of n x n 

symmetric positive-semidefinite matrices. 

{a'fi} M : kl X [0,T] ^ R"^"- is uniformly bounded and takes values in the space of n x n 
symmetric positive-definite matrices. 

(og) ^ : rZ —R”^” Zs bounded, -measurable and takes values in the space of n x n 
symmetric positive-semidefinite matrices. 

{a'rj) A*, 7 * : rZ X [0,T] ^ R for i G {1, • • • ,n} are uniformly bounded and strictly positive. 

(a'g) /3 : rZ X [0, T] —)■ is uniformly bounded and takes values in the space of n x n 

symmetric matrices. 

(a'g) There is no simultaneous jump among (A/"*, 

5.1 The market description 

We consider a market quite similar to what is described in Section 1,1.11 but now with n 
securities. The market maker’s position for the securities starting x G R” at time t is given 
by the following n-dimensional vector process: 

Tl ng n ng TL ng 

X'f’^{t,x) = X -\-'S^ / / eizJ\f^{du,dz) -\- / 7 r„du + ^ (5.1) 

^ Jt Jk Jt Jt 

where ej, i G {I,-- - ,n} is the unit R”' vector whose elements are all zero except the i- 
th element given by 1. tt = ... and 6 = (A)jgp ... are F-predictable trading 

strategies of the market maker in the exchange and in the dark pool, respectively. The 
superscript i is added to distinguish the corresponding security. Wfi G {I,-- - ,n} is the 
counting process denoting the occurrence of the execution of the Z-th security in the dark pool. 
AA*,z G {I,-- - ,n} is the counting measure which describes the occurrence of an incoming 
customer order of the i-th security and its size. 

We suppose that there exists a common bounded support K C R\{0} for the size of the 
incoming orders. We assume as before the existence of the compensators such that 
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for t G [0,T] are F-martingales for every i G {1, •'' Let us also set the stochastic 

processes $ = ^ and $2 = (^ 2 )ie{i,-,n} representing the 

moments of the order size by 

:= f zA'^{t, z)dz, 'I't •= [ \z\A'^{t,z)dz, [ z‘^A^{t,z)dz , (5.3) 

Jk Jk ' Jk 

for t G [0,T] which are all uniformly bounded by Assumption A'. 

We assume that the price vector S ’ {t,x) = (S'L‘^(t, which denotes the 

market price observed in the exchange under the impact of the market maker’s strategy 
(tt, S) starting from the position size x at time t, is given by 

Sf{t, x) = Ss + MsTTs - x) (5.4) 

for s G [t,T]. Here, M and /3 are not necessarily diagonal and hence they can induce direct 
as well as contagious stochastic linear price impacts from the continuous trading and also 
from the aggregate reactions of the other investors regarding the inventory size of the market 
maker. We can naturally imagine that, for example, due to the proxy hedging by correlated 
assets, a high trading speed or a big outstanding position of a certain security induces similar 
price actions among the closely related assets. 


5.2 The market maker’s problem 

We model the cash flow in the interval ]t,T] to the market maker with strategy ( 7 r,< 5 ) as 


ds-^ [ [ Sif_{t,x){l-sgn{z)bl)zJ\f\ds,dz) 

i=idi Jk 

Y, I USs- f3sX::f{t,x)ySl + 7rs\dT)dHl + r ljx:’yt,x)ds, (5.5) 


Sg’ {t, x) TVgds — 

n ^.rp 
t 


i=l 


where the symbol T denotes the transposition. We consider the following market maker’s 
problem: 


V(t,x) = ess inf E 

7r,Seu 


{xyy^xy^ + / {xyy^xyds 


+ [ {iSy^V'^s-lJxy)ds + Y [ [ Syy_{l-sgn{z)bl)zAf\ds,dz) 

Jt ^ Jt JK 

+E / ({Ss-^sxyy6i+?rg\5y)dHi 


i=l 


n 


(5.6) 


where we have omitted the argument (t, x) to save the space. By making ^ and 7 proportional 
to the (stochastic) covariance matrix among the securities, the market maker can include 
the portfolio diversification effects. In the above modeling, the customer orders and the 
executions in the dark pool are assumed to occur independently for each security. However, 
it is not difficult to introduce simultaneous customer orders or the dark pool executions 
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for an arbitrary subset of the securities by following the idea of dynamic Markov copula 
model studied by Bielecki, Cousin, Crepey &: Herbertsson (2014a, 2014b) [THlttl]. If there 
exist strong clusterings among the customer orders or the executions, an extension to this 
direction may become worthwhile. 

The set of admissible strategies U is defined below. 

Definition 5.1. We define the admissible strategies U by the set of ¥-predictable processes 
(tt, 5) that belong to IHI^(0, T) x ]H^(0, T) and also Markovian with respect to the position size, 
i.e., they are expressed with some measurable functions {f^,f^) by 

TT, = r{s,x:;^it,x)), S, = fis,x:;^it,x)) (5.7) 

where, for a € { 7 r,( 5 }, /“ : 11 x [0,T] x M” —)• M"" and /“(•,*) is an ¥^-adapted process for 
all a; G R”. 

Let us write the dynamics of the F'^-adapted bounded process fi as 


d 

djSt = + '^{a^).dWl 

i=i 


(5.8) 


and define 


7:=7 + y 


rf ■=rf' + 


(/3)m 


for f G {1, • • • , n} . 


(5.9) 


Assumption B' 

{b'^) pL^ : VL X [0,T] ^ and {(T^)i,i G {I,-- - ,d} : 11 x [0,T] —>• R"’^"' are uniformly 

bounded, symmetric and -adapted. 

{b' 2 ) C is positive-semidefinite. 

(63) 7 is positive-semidefinite dF 0 dt-a.e.. 

(64) There exists a constant c > 0 such that Xrf > c dF ® dt-a.e. for every 1 G { 1 , • • • ,n} 
and also My > c|yp dR ® dt-a.e. for every y G R"". 

Definition 5.2. The cost function for the market maker with a given position size x G R” 
at t & [0, T] is 


Jt-^(7r,(5) =E 


+ (A-’'5)^(/3,(6^), - ls))ds 


+ 


r I ttJM.tt, + {Ss + + Sl@s + jz (a: [r?:(d:)2 + + 

I i=l 




4 > 


2,s 


ds 




(5.10) 


where 0, b'if : H x [ 0 , T] —R” are defined by (0*)* := 4*5 — byi and {b^)l = b^l for 
z G {1, • • • ,n}. Here, the argument {t,x) of the position size is omitted to save the space. 
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Proposition 5.1. Under Assumptions A' and B', the market maker’s problem i5.6\) is equiv¬ 
alent to 

V{t,x)=ess inf (5-11) 

{n,S)GVl 

and it has a unique optimal solution gU. 

Proof. By using 

+ £ (^^iX:’Yu^sX:’^ds + . dWs 

Ti / 'fp 'fp \ 

+ E(/ +jJ^{fi,),,z^M\ds,dz)^ 

and redefining the value function 

V{t, x) := V{t, x) - ^x^Ptx (5-12) 

one can prove it in exactly the same way as Proposition 13.11 □ 


6 Solving the problem with multiple securities 


6.1 A candidate solution 


We derive a candidate solution for the market maker’s problem. Firstly, let us rewrite the 
optimality principle for the problem with multiple securities. 

Proposition 6.1. (Optimality Principle) Let Assumptions A' and B' are satisfied. Then, 
(a) For all x € M”, (7r,5) G U and t G [0, T], the process 


P(s, - Q)du + I' 

+ {Su + + S:,Si] + 


TV'). 


^l,u)\du\ ( 6 . 1 ) 

> / se[t,T] 


is an ¥-submartingale. 

(b) {n*,S*) is optimal if and only if 


P(s, Af A*) + J' ((Af’'5*)T7,Af’-5* + (Af’^‘)T - Q^du + | (O^MJO 


+ (S, + MAVK + Sl@u + f2{^u[vi{K? + S16*:,] + 


du 


( 6 . 2 ) 


i=l 


s(^[t,T] 


is an ¥-martingale. 
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Derivation of a candidate solution and the associated stochastic HJB equation is similar 


to the single security case. We assume that the semimartingale 
following decomposition: 


X 


iG[0,T] 


has the 


rs ns 

V{s,x) = V{t,x) + / a{u,x)du+ / Z{u,x)dWu 

Jt Jt 


(6.3) 


where a : D x [0, T] x M” —)■ R, Z : D x [0, T] x R” —)■ R*^ and a{-,x) as well as Z(-, x) are F^- 
adapted processes for all x G R"". We suppose that the value function can be decomposed, 
for every t G [0,T] and x G R” as 

V{t,x) = x'^V 2 {t)x + 2x'^Vi{t) + Vo{t) (6.4) 

Z(t, x) = x~^ Z 2 {t)x + 2x~^ Zi{t) + Zo(t) (6.5) 


where V 2 : Q x [0,T] ^ R’"''", Di : D x [0,T] ^ R’", Dq : x [0,r] ^ R, Z2 : D x [0,r] ^ 
]^nxnxd^ Zi : D X [0, T] —)■ and Zq : D x [0,T] —)■ R'^ are all F^-adapted processes. In 

addition, V2 and Z2 (with respect to the hrst two indexes) are symmetric. 

A lengthy but straightforward calculation shows that a necessary condition for the opti¬ 
mality principle is 


a{u, x) + x'^-fux + x'^{Pu {b^)u - In) + SI&U + Y1 

i=l 

n 

+2x'^V2{u)^u + 2Vi{u) 

i=l 


-hinfj (tv + M^^[V2{u)x + Vi{u)+ ^{Su +Mu&u)]J Mu 

X (tt + [V2{u)x + Di(«) + i(S,. + M„ 0 „)]) 

\V2{u)x + Vi{u) + ^Su]~^ 

[V2{u)]i,i + Tji 


+ X/ + du) (<^* + 


2=1 


V 2 {u)x + Vi{u) + -{Su + Mu&^ 


iT 


M, 


-1 


V 2 {u)x + 1/1(11) + i(s„ + M„0 


2=1 


([V2{u)x + Vi(tt) -I- - 

[V'2(«)]v + d\ 


e* 


> = 0 dF ® dt — a.e.. 


( 6 . 6 ) 


where we need [V2]i,i -|- 77* > 0 dP (8* dt-a.e. for every i G {1, • • • ,n}. As a result, we obtain 
the following. 
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A Candidate Solution 

A “candidate” of the optimal solution and the corresponding value function for the market 
maker’s problem i5.11\) are given by 

< = -M-^(v 2 iu)XlT {t,x) + V,{u) + ^{Su + MuSu)) (6.7) 




[v2{u)x:y\t,x)+v,{u) + -SuY 


[V2iu)]i,i + pi, 


for z G {1, • • • , n} (6.8) 


for u ^ [t,T] and V{t,x) = X V 2 {t)x + 2x Vi{t) + Vo{t), respectively. Here, X'" ’ {t,x) is 
the solution of 

n 

Xf'^\t,x) = x + ^ 


2=1 


t JK 


ps ^ ps 

izM\du,dz) + / 7r*(iM + ^ / ei{5l,ydHl„ sG[t, r](6.9) 

Jt Jt 


(V2, Z 2 ), iVi, Zi) and (Vq, Zq) must be the well-defined solutions of the following three BSDEs 


V2{t)=f + j^ {-V2{U) 


^ +diag(^ 


A,, 


V'2 (m) + Pu 

A 


V 2 (u) du- Z 2 {u)dWu 

( 6 . 10 ) 


yi{i) = -[ \ V 2 iu) ^+ diag(———) 

Jt I ^V2[u) + Pu'' 


+ V2iu) 


^ +diag(^ 


A. 


V2(u) + PuJ_ 

Su 1 


V2{u) + Pu 


Vi{u)--{(3u{b’^)u-lu) 

■T 


2 -^eu-{b^)uj'ldu- Zi{u)dWu 

( 6 . 11 ) 


Voit) = - 


u) + 


2 

,T. 


+ diag(^- 


u) + 


V2{u) + PuJ _ 

+ {b^)u) ' n{u) - 

2=1 

i(5;[0„ + ^[l3n]iM,u) + \&uMu&u \du- [ Zo{u)dWu (6.12) 


2=1 


satisfying, for every i G {1, • • • ,n}. 


[V2]i,i + ry' > 0 


(6.13) 


dF®dt-a.e. inklx [0,r]. In the above, diag ^j is defined as a diagonal matrix whose 

A* 

(z, i)-th element z G {1, • • • , rz} is given by 


[V2iu)]pi + p. 
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6.2 Verification 


In the multiple-security setup, V 2 follows a non-linear matrix valued BSDE. Since there is 
no comparison theorem known for a multi-dimensional BSDE in general, we cannot apply 
the technique used in the single-security case. Interestingly however, we shall see V 2 is the 
backward stochastic Riccati differential equation (BSRDE) associated with a special type of 
stochastic linear quadratic control (SLQC) problem in a diffusion setup studied by Bismut 
(1976) dg. 

Theorem 6.1. Under Assumptions A' and B', there exists a unique solution of{V 2 ,Z 2 ) for 
the BSDE 116.10\) . In particular, V 2 takes values in the space of n x n symmetric positive- 
semidefinite matrices and is a.s. uniformly bounded i.e., there exists a positive constant C 
such that 

ess sup [ sup |V 2 (i)|(w) I < C , (6-14) 

\te[o,T] J 

and Z 2 G P > 

Proof. Let us introduce an n-dimensional Brownian motion w which is orthogonal to W and 
consider Ej- := V Ef’ where Ef^ is the augmented filtration generated by w. We study an 
n-dimensional (F' := -adapted vector process staring from x G M” at time t which is 

controlled by the 2n-dimensional vector process 6: 

X^^{t,x) = x+ [ CuOudup'^ [ DlOudwl , s e [t,T] . (6.15) 

Jt Jt 


Here, C : D x [0, T] —)• jg (defined by 

Cu ■= (inxn diag(A(,)) (6.16) 

for u G [0,T], where Inxn is the n-dimensional identity matrix, and diag(A*) is the n- 
dimensional diagonal matrix whose {i, i)-th element i G {1, • • • , n} is given by A*. We use the 
same notation for the diagonal matrices below. D* : D x [0,T] —?■ ^ ^ ... 

has zero entry for all except the (i, n -|- i)-th element which is given by 


[Dl 


u\i,n+i 



(6.17) 


for u G [0, r]. We define the admissible strategies U' as the set of 2n-dimensional F'-adapted 
processes 6 that belong to Bl 2 „( 0 , T). 


Now, let us consider the following SLQC problem: 

rT 


V'it, x) = ess inf E 

e&w 




Ei 


(6.18) 


where the argument (t, x) is omitted from X to save the space, and N : D x [0, T] 


p2nx2n 
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is defined for u E [0, T] by 


N„ = 


Mu 0-7 y 

Onxn diag(A^r?J,) J 

Then, by Proposition 5.1 in [12], the associated BSRDE is given by 

T ( ^ T 

P{t)=^ + j l-Piu)Cu[Nu + ^{DiyP{u)Diy^CjPiu)+jAdu-J Zp{u)dWu 

( 6 . 20 ) 

where P is connected to the value function as V'{t,x) = x~^P{t)x. Note that the stochastic 
integration by dw vanishes because W P w and that the terminal value ^ and all the processes 
included in the driver are F'^-adapted. By noticing that 


Nu + P{u)Dl 


i=l 


Mu Onxn 

Onxn diag(A^([P(u)]i,i + r]\ 


( 6 . 21 ) 


one can confirm that the BSDE of P is equal to that of V2 given by ([R 


Under Assumptions A' and B', ^ is positive-semidefinite and bounded, 7 is positive- 
semidefinite and uniformly bounded, C, D and N are uniformly bounded. In particular, 
there exists a constant c > 0 such that 


NuU > c\y\^^ dF®dt — a.e. 


( 6 . 22 ) 


for all y E Thus, by Theorem 6.1 in [12j . P (and hence V2) has a unique solution, which 
is symmetric, positive-semidefinite and a.s. uniformly bounded. In particular, this implies 
[P{u)\i^i >0, dP (8) dt-a.e.. 

Since P is positive, one sees from (I6.2ip . 


0 < (Nu + Y,iDiVPiu)Di) \ < 


2 = 1 


\y\ 


dF ® dt — a.e. 


(6.23) 


for all y E and hence (Nu + Yl'i=i{Pu)^P{'P)Pli\ is a uniformly bounded linear 

V /ug[0,T] 

operator. Using the boundedness of P and the other variables, one sees 


mt = P{t) - P{0) - j \ P{u)Cu(Nu + AZp{u) -^u\du (6.24) 


2 = 1 


for t E [0,T] is a uniformly bounded martingale. Thus, from the BDG inequality, for any 
p > 0, there exists a positive constant C such that 


E 




0 


< CE 


m 


< 00 


(6.25) 


°It is not difficult to confirm the same BSRDE arises as the stochastic HJB equation by the same method 
we used. 
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and hence Zp (and so does Z^) belongs to '^'P ^ ^■ 


□ 


For more general resnlts on the SLQC problem and the associated BSRDE, we refer to 
Peng (1992) [38] and Tang (2003, 2014) jlSlIll], where the assumption of the orthogonality 
“u; T W” is removed. 

The following results are obtained in exactly the same way in Proposition 14.3114.41 and 
Corollary 14.11 


Proposition 6.2. Under Assumptions A' and B', there exist unique solutions (Vi,Zi) E 
S^(0, T) X EI^^^(0,T) for h6.11\) . and {Vo,Zo) E S^(0,T) x 11^(0, T) for i6.1S\) . respectively. 

Furthermore, the process for the position size (X 




sG[t,T] 


given by 116. belongs to 

S^(t,T). The candidate solution ( 7 r*,< 5 *) given by and 116.8\) is well-defined and satisfies 

eSi{t,T) xSUt,T) CU. 


The above results establish the main theorem. 


Theorem 6.2. Under Assumptions A' and B', the candidate solution {it*, 6*) given by B) 
and i6.8\) is, in fact, the unique optimal solution of the market maker’s problem given by 

(KW- 


Proof. The proof is the same as that of Theorem 14.11 


□ 


Remark: A determination of the bid/offer spreads 

Before closing the section, let us comment on a possible determination of the bid/offer spreads 
b. Although we have assumed that the market maker do not dynamically control the bid/offer 
spreads to give a bias to the order flows, it is important of course to use a sustainable spread 
size for its market making business. Suppose, for example, the spread size 5* is proportional 
to the volatility |(T*| of the z-th security as 

K = aksi (6.26) 

where a > 0 is some constant and z E {1, • • • , n}. Even if the intensity (and/or distribution) 
of the customer orders is a non-linear function of ( 6 *)jg{i^... the market maker can obtain 
the cost function or the distribution of its revenue by running the simulation based on the 
optimal strategy ( 7 r*,( 5 *) for each choice of a, which will give enough information to fix the 
size of d. 


7 Implementation for a simple case 

In this section, we discuss the evaluation scheme for a simple case where V 2 becomes non- 
random. As we shall see below, the implementation of the optimal strategy is quite simple 
in this case. 

Consider a setup where 7 , M, 7 , A (and hence naturally so is /3) are non-random. In this 
case, V 2 is a solution of the following matrix-valued ordinary differential equation (ODE): 

^ = U(s) (m-‘ + diag(^^A_)) V 2 (s) -I., s€ \t,T] (7.1) 

V2(r) = { (7.2) 
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which is the same ODE studied by Kratz &: Schoneborn m- As long as 7 , M, 7 , A satisfy 
the boundedness conditions in Assumptions A' and B', this Riccati equation has a positive 
bounded solution. It is not difficult to numerically solve this equation by the standard 
technique for ODEs. In contrast to the model in m , we still need to evaluate Vi to implement 
the optimal strategy (See Eqs (|6.7I) and (|5.8p .l. 

Eor notational simplicity, let us put 

F{s) := V 2 is) + diag(^^r^-^)) , s E [t,T], (7.3) 

which is a deterministic matrix process. Let us also consider another deterministic matrix 
process Yt^. defined by the ODE 


dYt,s 

ds 


F{s)Yt,s, sG[t,T] 


where Yt^t = Inxn- Then, we have 


ds 


-Y,-^F{s), se[t,T] 


(7.4) 


(7.5) 


with Y^j- = Inxn and it is straightforward to obtain 


V,{s) =-Yt,s 


^F{u)Su-V2{u)(^^®u + {bAf\ 


- -{PuibAf)^ - lu) 


F 


w 


du 


(7.6) 


for s E [t,T]. One sees that the bid/offer spread, the customer order flows and the size of 
repo rate impact the optimal strategy though Vi. Its evaluation only requires 


E[$„|jf], E[( 6 ^')JJ-f], E[lu\F 


wi 


(7.7) 


These quantities can be obtained analytically for simple models. Otherwise, one can apply 
the standard small-diffusion asymptotic expansion technique, which is developed by Yoshida 
(1992a) [l 6 ], Takahashi (1999) [42], Kunitomo & Takahashi (2003) [32| for the pricing of 
European contingent claims, and also Yoshida (1992b) [IT] for statistical applications. See 
Takahashi (2015) [48] and references therein for the recent developments. 


8 Implementation for a general Markovian case 

Although it is impossible to solve V 2 analytically in a general setup, getting an explicit ex¬ 
pression of its approximation is very important for successful implementation of the proposed 
scheme. A similar BSDE is also relevant for solving a different type of optimal liquidation 
problem treated in [5]. Eurthermore, considering the wide spread applications of SLQC prob¬ 
lems in various engineering issues, developing a successful approximation scheme for a general 
BSRDE should be a very important research topic in its own light. 

There exists an analytical approximation technique for non-linear BSDEs, which was 
proposed in Eujii & Takahashi (2012a) [22]. The method introduces a perturbation to the 
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driver and then linearizing the BSDE in each approximation order. It then adopts the small- 
diffusion asymptotic expansion to evaluate the resultant linear BSDEs. Its justification for 
the Lipschitz driver was recently given by Takahashi & Yamada (2013) [IS]. We refer to 
Fujii &: Takahashi (2012b) [23| for an example of explicit calculation, and Fujii & Takahashi 
(2014) [2T| as efficient Monte Carlo implementation where the analytical calculation is too 
cumbersome. See also Shiraya & Takahashi (2014) [IT] and Crepey &: Song (2014) [T7| as 
concrete applications of the proposed perturbation method to the so-called credit valuation 
adjustment (CVA). 

In this section, we propose a different type of perturbative expansion method. In contrast 
to the method [ 22 ], it does not require the perturbation of the driver and allows simpler 
analysis. The new method directly expands the BSDE around the small-diffusion limit of 
the associated forward SDE, which only yields a system of linear ODEs to be solved at each 
order of expansion. 


8.1 A perturbative expansion scheme 

Let us first explain the idea of our perturbation scheme. We will provide the justification and 
error estimate later. For clarity of demonstration, let us assume a single security case. The 
extension to a multiple security case is straightforward. We introduce the underlying factor 
process A ; D x [t, T] —>■ with an arbitrary starting time t G [0, T], which follows the SDE 
(do not confuse it with the position size process): 


A 


s 


= X + 


^{u, Xlf)du + 


a{u,X^^^)dWu . 


( 8 . 1 ) 


where fi : [t,T] x —>■ and a : [t,T] x —)• The superscript {t,x) indicates the 

initial condition for the process, which will be omitted if it is clear from the context. Let a 
function / : R'^ x R —R be defined by 


f{x,v) := - 


1 




A(x) 


M{x) V + ri{x) 
where M, 7 , r/, A : R"^ —>■ R, and consider the BSDE 


-t- 7(x) 


Jt Jt 


( 8 . 2 ) 


(8.3) 


where V : D x [t, T] —>■ R, Z : D x [t, T] ^ R'^^'^. This BSDE corresponds to (|4.11l) with 
stochastic coefficients driven by a Markovian factor process A. 


Assumption P 

1. The coefficients fi,a are bounded Borel functions, and ffit,x) and a{t,x) are continuous 
in (t,x) and smooth in x with bounded derivatives of all orders. 

2. There exist constants 01,02 > 0 such that for My G R'^ andM{t,x) G [0, T] x R'^, 

ai|y|^ < y^[aa'^]{t,x)y < 02 |?/|^. (8.4) 


29 






3. ^{x), M{x),'y{x),ri{x), X{x) are bounded smooth functions o/x G and satisfy Assump¬ 
tions A and B with Vx G They are also assumed to have bounded derivatives of all 

orders. 

Now, in order to approximate this process, we introduce a small parameter e g]0, 1] and 
e-dependent process X'^ 

= p.{s, Xl'^'^)ds + €a{s, Xl’^’^)dWs s G [t, T], xj’^’" = x . (8.5) 


By using this process, we carry out small-difFusion asymptotic expansion of the system. The 
associated perturbed BSDE is given by 



At the moment, let us assume the differentiability in terms of e so that we have 


= Xpl + cAW + e^Apl + ... 

y-t,x,e ^ -^jo] ^ ^yjl] ^ ^ . 

^t,x,e ^ ^[0] ^ ^^[1] ^ ^2^[2] ^ . 


for s G [t,T], where we have defined 


yH . 

^ S 


1 

n! (9e^ 



e=0 


s G [t, T] 


and similarly for the others. 

For the zero-th order, and are given by the solutions of the ODEs; 


(8.7) 


( 8 . 8 ) 


= sG[t,T], Af=x 

= -/(Apl,Fio]) se[t,T], fM=^(40]) 


(8.9) 

( 8 . 10 ) 


which corresponds to a deterministic case discussed in the previous section. Thanks to 
Assumption P, the above Riccati ODE has a positive bounded solution. We obviously have 

ZM = 0. 

In the first order, we have the linear FBSDE system: 

dAW = d^n^{s)xl^Us + a^{s)dWs s G [t,T], A^ = 0 (8.11) 

fJ'I = a,^0(r)AW + I 5„/°(s)Fiil + 5 ,/°(s)aW ^ ds - 

( 8 . 12 ) 


where we have used the short hand notation; 

/(s) := m(a 4’), := <^(^,4°’), C°(r) := ^xP), f\.s) := /(Af ,4') • (8.13) 

which are all deterministic functions. We have also used dx '■= {d/dxi)i<i<d, := d/dv. 
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By the assumptions we made, it is clear that there exists a unique solution for the BSDE 
satisfying G 8>P{t,T), G MP(t,T) for Vp > 1. In fact, it is straightforward to explicitly 
solve it as 


EW=y(s)^xW, sG[t,T] 


where y : [t, T] is the solution of the linear ODE: 


(8.14) 


J=1 

y{T) = d,e{T). (8.15) 


Due to the assumption on boundedness, = y{s)^ {s) is actually uniformly bounded. 

In the second order expansion, one can find 

dxf^ = (dy{s)xf^ + i(xW)^(92/(5))xW)ds + xW5,a°(s)dIE„ s G %T] 

Xf = 0 (8.16) 


= d^e{T)x'^? + \{X)^V {die{T))X^^^ + j^^{d,fis)VP + 9,/°(s)xl2] 

+ld^J\s){vmf + d,,,f{s)xl^wm + ^{xl^Y{dlf%s))xm}ds - [ zf^dWs, 

(8.17) 


which is also a linear BSDE. In this case, one has the solution of the following from: 

1"?’ = y 2 {s)^xf^ + (xW)Tyi(s)xW + y^{s), s G [t,T]. (8.18) 

Here, y 2 : [t,T] —)• M'^, yi : [t,T] and yo : [t,T] —)■ M are defined as the solution of the 

next linear ODE system for s G [t,T], i,j G {1, • • • ,d}: 


d[y2{s)]i 

ds 


■(^5xi(/(s))i[y2(s)]i + dyfis)[y2{s)]i'^ - d^,f{s) 


i=i 

d 


+ [yiis)]jAif^\s)%] +dyf{s)[y,is)],,, 

^ k=l 


'^[y 2 {s)]kdxi,xj{y°(.s))k - ^d‘^f{s)[y{s)]i[y{s)]j 
^ k=l 

-\{[y{s)]jdxi + [y{s)]idx^dyf{s) - ]^dxi,xjf{s) 

= -dvf{s)yQ{s) - Tr[yi(s)fT°(s)(o-°(s))’^] (8.19) 
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with the terminal conditions: 

y 2 {T) = da,e{T), [yi{T)]ij = ^d.„..e{T), yo{T) = 0. (8.20) 

The above results can be obtained by applying Ito formula to (I8.18p and comparing its drift 
to that of (I8.17p . See Fujii (2015) |2T] as a related idea of expansion of BSDEs. 

These procedures can be repeated to an arbitrary higher order. In each order, one can 
show that there exists a unique solution E §P{t,T), Zl”! E EI^'(t, T) for Vp > 1 due to 
the boundedness assumptions and the linearity of the BSDE (See Theorem 5.17 in [37].). 
The solution can be expressed by some polynomials of {-TW}i<j<„ with the coefficients to 
be determined by the linear ODE system. Note that, for every order n, E §P{t,T) for 
Vp > 1. 


8.2 Convergence 

Theorem 8.1. Under Assumptions A, B and P, there exists some positive constants C,C', 
which are independent of e, such that 


E 




N 

^y[o] 

n=l 


P 

%T] 




( 8 . 21 ) 


( 8 . 22 ) 


for Vp > 1 and every positive integer N. 

Proof. The arguments for the justification and convergence are similar to those of |45j and 
we sketch them in the following 0. Under Assumption P, the continuity and differentiability 
of with respect to e are well known (See, for example, jUj.). For the continuity and 

differentiability of we can follow the same arguments of Section 2.4 of El Karoui 

et.al. (1997) [T9| and Theorem 3.1 of Ma &: Zhang (2002) |3l|. Their results are based on the 
popular estimate (See Lemma 2.2 of [34]) for a BSDE with a Lipschitz driver. Although the 
driver / is not Lipschitz in our case, we can in fact use the same estimate. This is because, 
we know the solution is uniformly bounded and non-negative thanks to Proposition 14.21 
For higher derivatives (cl^E^,(9^Z^)fc>i, the arguments are more straightforward since the 
BSDEs are linear for them. One can check that the Lipschitz conditions are satisfied by the 
similar reasons. Thus, recursively, the arguments in [laiM] guarantee the continuity and 
differentiability up to an arbitrary order. 

Now consider the n-th order derivatives 


an an 

A _wi.a:,e A_ 

n,s ■ Q^n s i ^n,s ■ Q^n s ) 

and their restriction to the condition e = 0 


T/N = 
s n' 


e=0 


ZN = Iz*;"/ 


n! 


e=0 


ytyX^e 

■^n,s 

f)n 

— Q^n s 

(8.23) 

xiA 

^ yt,x,e 

- nr^ e=o' 

(8.24) 


^The author is grateful to prof. Takahashi for helpful discussions. 
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Taylor expansion gives the associated BSDE as (omitting superscript (t,x)), for s G [t,T], 


— Gn + 


Hn,r + dr - Z^.dWr (8.25) 


where 


n ^ 

G» ="! E E n i3iX,.T 

k=l /3iH-l-/3fc='n,,/3i>l j=l ^ 


(8.26) 


and similarly 


k—i 


^n,r ^ 


!E E E E 

k=2 /3iH ^=0 j=k—i-\-l 

7 = 1 ^ 1 1- ^ , 1 




i\{k — i)\ 


(8.27) 


One obtains the SDE for X^ ^ in a similar manner. For every n, due to Assumption P 
and the linearity of the SDE, one can show that X^ . G §^(t, T) for every Vp > 1 and e g]0, 1]. 

Thus < oo holds for Vp > 1. Then, since (j8.25j) is a linear BSDE, one can show 

recursively that \Hn,s\ds^ < oo for Vp > 1 and that there exists a unique solution 

G §P{t,T), . G 'EP{t,T) for every n, Vp > 1 and Ve g]0, 1], by applying Theorem 5.17 


in [37] 

Using Taylor formula, one sees 


N 


_ ,n o? 


n=l 


n! de'' 


E=0 


+ e 


.AT+l 


I 


^ (1 - u)^ (9^+^ 


= Ujol + ^ e-uN + ^ / (1 - 

n=l do 

Thus, there exists some constant C such that 


N\ 

du. 


VI 


du 


(8.28) 


N 


E 

and similarly. 


V^- 




n=l 


[i.T] 


< ,pX+i)C VW+I, I \l^T] (8-29) 


E 


^1- 


+ 


N \ pG 

^ds] 


< e 


,p(iv+i)(^E ( 


■T . rl 


n=l 




\Jt Oo 


\ZW+i,s\‘^du\ ds] 


N p/2 

J 


□ 


for every N and Vp > 1, which proves the claim. 

The above result can easily be extended to the multiple security case by using Assumptions 
A', B' and the boundedness of V 2 proved in Theorem 16.11 Once the terminal penalty is 
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replaced by a random variable ^ G LP(n), the proposed perturbation algorithm can be applied 
to a different class of BSDEs [5], too. Detailed numerical tests and the extension to more 
general class of BSRDEs [13] will be left for an important future work. 

9 Concluding Remarks 

In this paper, we discussed the optimal position management strategy for a maker maker who 
faces uncertain in- and out-flow of customer orders. The optimal strategy is represented by the 
solution of the stochastic Hamilton-Jacobi-Bellman equation which is decomposed into three 
(one non-linear and two linear) BSDEs. We provided the verification of the solution using 
the standard BSDE techniques for the single-security case and an interesting connection to a 
special type of SLQC problem for the multiple-security case. We also proposed a perturbative 
approximation technique for the relevant BSRDE, which only requires a system of linear 
ODEs to be solved at each order of expansion. Its justification and error estimate were also 
given. 

Assuming general F-adaptedness (instead of F’^-adaptedness) of the relevant parameters 
looks an interesting extension of the proposed framework. This situation will arise when one 
introduces simultaneous jumps in the parameters, such as M, and the executions in the dark 
pool. In this case, the driver of the resultant BSRDE depends on the martingale coefficient 
of the counting process. As long as we know, the existence and uniqueness of the solution for 
the corresponding BSRDE have not yet been proved. 

It looks also interesting to combine a stochastic filtering for the intensity of customer 
orders. Introducing a hidden Markov process, for example, is likely to help to model possible 
herding behavior among the customer orders. See a related work Eujii Sz Takahashi (2015) [25| 
on the mean-variance hedging problem for fund and insnrance managers. 
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A The Property of the Terminal Position Size 

In this appendix, we study the behavior of the remaining inventory Xj, ’ at the terminal 
time according to the change of the penalty size. We shall prove that it can be made arbitrary 
small by increasing the size of the penalty This result implies that the proposed strategy can 
be considered as a generalization of the optimal liquidation solution in the existing literature 
to the situation with uncertain customer orders. 

Let us take a positive constant 1 < L < oo and set ^ = L, i.e., ^ = (5t/‘^ + L. We denote 
the corresponding solutions of the BSDEs (j4.1ip . (I4.12p and (I4.13|) by (Ej^,-^/'){j=i, 2 , 3 }) 
respectively. 
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Assumption C 

Take the lower bound c in the assumption ( 63 ) in sueh a way that c/(l + A) < 1 and also 
c := c/[M(l + A)] < 1/2. Obviously, one can always choose c > 0 (or equivalently M,X) to 
satisfy these inequalities. 


Lemma A.l. Under Assumptions A, B, C and = L, the following inequalities hold for 
every 0 < t < s < T with an L-independent positive constant C; 


V^t) < C 
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T — t -\- Cl 


exp 




:n < 


- r{u,V 2 {u))d 

andr{t,y) 


T - S +CL 
T — t + Cl 


(A.l) 


i y + vt 


Proof. Since the inequality in Proposition 14.21 holds arbitrary e > 0, one can choose e = cl = 
1/L. Then one obtains 
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(A.2) 


Similarly, 
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(A.3) 


where c > 0 is the lower bound given in ( 63 ). Thus, 
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It yields 
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Note that for every 0 < t < s < T, 


( xeL+T—s 
y xe^+T—t 
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is a increasing function for x > 0. Thus, due 
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to the arrangement of c, one obtains 


exp 



r{u,V^{u))du^ 


/ T - s + ei 
\T — t + eL 


C 


(A.6) 

□ 


We also have the following Lemma. 

Lemma A.2. Under Assumptions A, B, C and ^ = L, there exists an L-independent positive 
constant C such that 
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Proof. Let us put A : n x [0, T] —>• M and a : 0 x [0, T] ^ M as 
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Obviously, A, a G S^(0,T) C S^(0,T), and whose S^(0, T)-norms can be dominated by L- 
independent constants. It is straightforward to check that can be written as 


= -E 
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Thus, by Lemma lA.ll it satisfies the following inequality for Vt G [0,r]: 
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Notice that (mt := E ||a||T + ||^||t|-T/^ ) is a square integrable martingale. Thus, 

V 1 J/te[o,T] 

from Doob’s maximum inequality, one has 



< CE 

1 |2 
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where the right-hand side can be dominated by an L-independent constant. 
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Now, let us define another process G : Q x [0, T] —R as 




V2 ('w) + h« 


5, 




0 . 


(A.13) 


which satisfies G € S^(0,T) C S^(0,T) and its S^(0,T) norm can be dominated by an L- 
independent constant by the first part of the proof. From (|4.8I) . (14.9p and (I4.10p . it is easy 
to see that 


X* =e-^0<^y2G))dUx_ / 


+ r / (A.14) 

iO JK Jo 


holds for every t E [0,T] (We used the notation X* := (0,x).). Using the fact that 

r(-, V^{-)) is a positive process and the BDG inequality, we have, with some L-independent 
constant G, 
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Let denote an L-independent constant dominating the first term by G'. Since we already 
know X* E §^(0,r). 
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and hence by the Gronwall lemma, 
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Combining the first part, the claims were proved. 


□ 


Then, we can establish the following result. 


Theorem A.l. Under Assumptions A, B, C and ^ = L, there exists an L-independent 
positive constant C satisfying 
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T Cl 
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and hence one can make the terminal position size arbitrarily small by taking a large L < 00 
as the penalty. 
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Proof. From (jA.14h and Lemma lA.ll we have 
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Notice that the expectation in the second term is dominated by an L-independent constant 
by Lemma lA.21 Using the assumption 2c < 1, we have 


E 
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CL \ 


2c 


e^/ L-2C V i +CL/ 1 - 2c 

<C(^j (A.20) 

with some L-independent positive constant, and hence obtained the desired result. □ 


Remark 

Although we can discuss the limit of the singular terminal condition L —>■ oo as presented in 
[5], we can only apply their results to V^. For there appears a singular drift term which 
is expected to create a discontinuity at the terminal point. This makes the detailed analysis 
difficult to carry out. However, as the previous result shows, we can make the terminal 
position size arbitrarily small by selecting large enough L < oo as the penalty. Therefore, 
the proposed strategy can also be used as an effective liquidation strategy in the presence of 
incoming customer orders for a market maker. 

Although it is natural, even in a multiple-security setup, to imagine that one can make the 
terminal position size arbitrarily small by increasing the size of the eigenvalues of Although 
it is intuitively clear, it is difficult to prove since we do not have an explicit expression for 
the upper/lower bound of V 2 any more. 

Let us suppose, in the interval [T — e, T] with some constant e > 0, that M, 7, fi can 
be diagonalized by the common constant orthogonal matrix O. In addition, suppose the 
market maker stops accepting the customer orders and stops using the dark pool. Then, by 
considering the securities in the base 0^5 and the corresponding positions X, the market 
maker’s problem can be decomposed into n single security liquidation problems. In this case, 
V 2 := 0^V20 becomes diagonal process in [T — e,T] and Ui := interacts with the only 

one corresponding element of V 2 . In this special situation, it is clear that the position can be 
made arbitrary small by the corresponding optimal strategy thanks to the arguments made 
in the single security case. 
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